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Abstract 

Unlike classical bifurcations, border-collision bifurcations occur when, for example, a fixed 
point of a continuous, piecewise C 1 map crosses a boundary in state space. Although classical 
bifurcations have been much studied, border-collision bifurcations are not well understood. This 
paper considers a particular class of border-collision bifurcations, i.e., border-collision period- 
doubling bifurcations. We apply a subharmonic perturbation to the bifurcation parameter, 
which is also known as alternate pacing, and we investigate the response under such pacing 
near the original bifurcation point. The resulting behavior is characterized quantitatively by a 
gain, which is the ratio of the response amplitude to the applied perturbation amplitude. The 
gain in a border-collision period-doubling bifurcation has a qualitatively different dependence 
on parameters from that of a classical period-doubling bifurcation. Perhaps surprisingly, the 
differences are more readily apparent if the gain is plotted vs. the perturbation amplitude (with 
the bifurcation parameter fixed) than if plotted vs. the bifurcation parameter (with the pertur- 
bation amplitude fixed). When this observation is exploited, the gain under alternate pacing 
provides a useful experimental tool to identify a border-collision period-doubling bifurcation. 



1 Introduction 

In contrast with the usual C 1 context for classical bifurcations jSj, the context of border-collision 
(B/C) bifurcations is continuous, piecewise C 1 maps (UEU]. srm pl es t normal form of a B/C 
bifurcation, which we derive below, occurs for iteration of a piecewise linear map as follows 

j , i = | Ax n + cfi, ifac^>0 ,^ 
\ Bx n + ciJ,, if Xn < 

where x is an m-dimensional vector, A and B are m x m constant matrices, c is a m-dimensional 
constant vector, and the scalar \x represents a bifurcation parameter. Here, the term border refers 
to x^ = 0, where the superscript indicates the first component of x: i.e., the border is the surface 
across which derivatives of the mapping jump. 

Examples of B/C bifurcations have been observed in various engineering and biological systems 
1131 1141 115j . Different phenomena in B/C bifurcations have been classified by a few authors 
PJ El HI 03 EE H| • Border-collision bifurcations may exhibit surprising phenomena such as so- 
called instant chaos, which is a direct transition from a fixed point to a chaotic attr actor |14j . 
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Figure 1: Schematic bifurcation diagram of a border-collision period-doubling bifurcation. Here, 
the thin horizontal line represents the border. Solid lines represent stable solutions whereas the 
dashed line represents the unstable fixed point. 

Alternatively, B/C bifurcations may bear great similarity to their smooth counterparts, such as 
B/C period-doubling bifurcations [H I2L)|. which are the main focus of this paper. 

Unlike a classical period-doubling bifurcation that occurs when one eigenvalue crosses the unit 
circle through —1, eigenvalues are not predictive for the onset of a B/C period-doubling bifurcation. 
Instead, as indicated in Figure ^ a B/C period-doubling bifurcation occurs when a branch of fixed 
point collides with a border. Despite such intrinsic difference, it may be hard to distinguish between 
the two bifurcations, especially based on necessarily limited data collected from experiments [2j|. 
This work investigates a possible technique to differentiate between classical and B/C period- 
doubling bifurcations, inspired by a phenomenon known as prebifurcation amplification previously 
observed in classical period-doubling bifurcations |17U18| . The term "prebifurcation amplification" 
refers to the fact that, near the onset of a classical period-doubling bifurcation, a subharmonic 
perturbation to the bifurcation parameter leads to amplified disturbances in the system response. 
In the literature, such subharmonic perturbations are also referred to as alternate pacing |2l ll2[ IT9]. 
The amplification due to alternate pacing of a classical bifurcation has been characterized by a gain 
relating the disturbances in response to the perturbation amplitude I n this paper, we 

compute the gain for a B/C period-doubling bifurcation under alternate pacing. Specifically, we find 
that the gain of a B/C bifurcation, as a function of the bifurcation parameter and the perturbation 
amplitude, differs significantly from that of a classical bifurcation. 

In Section 2, we first investigate the conditions for the existence and stability of a B/C period- 
doubling bifurcation. Then, the response of a B/C period-doubling bifurcation to alternate pacing 
is derived in Section 3. Section 4 quantifies the effect of prebifurcation amplification through a 
gain. Finally, Section 5 provides a concluding discussion. 

2 Hypotheses and Notation 

2.1 Normal form of a border-collision bifurcation 

We first derive the normal form (^Q). Consider iteration under a general border-collision (B/C) map 

z n +l = f {z n , v) , (2) 
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where z n G M m , v is a bifurcation parameter, and 

f / 7 ,^_ / /a^.^). if^O,z/)>0 

/( ' j "l /„(*,»), if/3(z,^)<0 • ^ 

To have a B/C bifurcation, we impose the following three conditions on the above map. 
Condition 1: Continuity on the boundary: f A (z, v) = f B (z, v) whenever (3 (z, v) = 0. 
Condition 2: For v = i/yf, © has a fixed point z\y$ on the boundary, i.e., 

P (zui, Vbii) = (4) 

and z hi f = fx {z\au ^bif) = /b Obif, ^bif ) ■ 

Condition 3: The boundary is locally nonsingular; in symbols 

V/3 | b if = V/3 (z bif , i/ bif ) ^0, (5) 

where V indicates the vector of derivatives of (3 with respective to zx, . . . ,z m ; i.e., it does not include 
the v derivative. 

To derive the normal form, choose an invertible matrix R s.t. 



rp 

where e\ = (1, 0, . . . , 0) . Let b be the vector 



R (V/3| bif ) =ei, (6) 
ictor 

b= (dvP\ hi{ ) e x . (7) 



Making a change of coordinates 

x = R (z- Zuf) + {v - fbif) b, ii = v- z^bif (8) 
and dropping all higher order terms in (J2J) yields a map that assumes the form Q, where 

A = R (Df A \ hi{ ) R-\ B = R(Df B \ m ) R- 1 (9) 
c = b-Ab-R (d„/ A | bif ) = b-Bb-R (d v f B \ hii ). (10) 

Here, -D/ denotes the differential of /, i.e., the mxm matrix of partial derivative of / with respective 
to zx, ■ ■ ■ z m . Note that the two expressions of the vector c are equal due to the continuity of /. As 
will be seen later, condition 4 below requires that c / 0. 

Remark: The continuity condition (#1 above) implies that Ax = Bx if x^ 1 ' = 0. Thus, only 
nonzero column of B — A is the first column. 



2.2 Conditions for the period-doubling case 

For the rest of the paper, we assume map JQ) undergoes a B/C period-doubling bifurcation as 
illustrated in Figure ^ Specifically, we assume that: 

Hypothesis 1: (a) For fi > 0, map (^Q) has a stable fixed point Xf p (/i) in the half space 
{a^ 1 ) > 0}. (b) For \i < 0, there exists a fixed point in the half space {xW < 0} but it is unstable. 

Hypothesis 2: For fi < 0, (^Q) has a stable period-two orbit, with one point in {a^ 1 ) > 0} and 
the other in ja^ 1 ) < 0}. 

Conditions 4-6 below guarantee that (^Q) possesses this structure. We do not assume that the 
fixed point for /i > is the only attractor nor do we assume the period-two orbit for fi < is the 
only attractor. Indeed, the analysis of Banerjee and Grebogi suggests that, more often than 
not, there will be additional attractors. 
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2.2.1 Analysis of fixed point solutions 

When fi > 0, a fixed point of © in the half space {x^ > 0} must satisfy 

X{ P (^) = ^X {p , (11) 

where Xf p is the constant vector 

X ip = (I - Ay 1 c. (12) 

By hypothesis 1 (a), we have the following condition. 

Condition 4: The first component of the vector Xf p is positive, i.e., X$ > 0, and all eigen- 
values of A are inside the unit circle, i.e., Vi, |Aj (A)\ < 1. 

When /x < 0, the unstable fixed point in {x^ < 0} is given by 

^unstable (A 1 ) = H{I ~ B)' 1 C, (13) 

and by hypothesis 1 (b) we have the following condition. 

Condition 5: The first component of the vector (/ — B)~ c is positive and at least one 
eigenvalue of B satisfies that \\ (B)\ > 1. 



2.2.2 Analysis of period-two solutions 

When fj, < 0, let us write x upP er (/u), Slower (m) f° r the stable period- two orbit of where as in 
Figure n ^upper (m) is the point in {x^ > 0}. These points satisfy 

^upper Slower ^upper (m) + /xc (14) 

It follows from the above equations that 

^upper — A^-^upperj Slower — A'-^lower) (1^) 

where 

^upper^-S^)" 1 (/ + £?) C, -Slower = (I — A B)^ 1 (I + A) c. (16) 

For these calculations to be consistent and for the orbit to be stable, we need the following condition. 

Condition 6: X$ pei < 0, X^ eT > 0, and Vi, |Aj(AS)| < 1. (Note that matrices AB and 
5^4 have same characteristic polynomials and hence have the same eigenvalues.) 



3 The response to alternate pacing 

Alternate pacing imposes a subharmonic perturbation to the bifurcation parameter [/,, rendering 
the map as 

= f Ax n + [ii + {-l) n 5] c, ifx^>0 
Xn+1 \ Bx n + [fi + (-l) n 5} c, ifx^<0 ' 

For use below, we define 

d = (I + A)' 1 c; (18) 
to avoid a degenerate response to alternate pacing, we impose the following condition. 
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Figure 2: Unilateral and bilateral solutions of alternate pacing (dashed) compared to the solutions 
when 5 = (solid). The thin thin horizontal line represents the border. 



Condition 7: The first component of the vector d is nonzero, i.e., 

We are interested in period- two solutions of l|17|). whose two points are denoted by y U pper (a 4 , S) 
and yiower (m, $)■ As depicted in Figure El and will be derived in the following, the character of 
the period-two solution depends on fi. 1.) For \x > /U cr i t (5), where A«crit (5) is defined in l|23j) 
below, the response is a unilateral solution: i.e., both points are above the border, or in symbols 
2/upper > yfowcr > 0> an< ^ ^.) ^or A 4 < Merit (<5)> the response is a bilateral solution: i.e., one point is 
above and the other is below the border, or in symbols j/upper > > yj^L_. 

3.1 The unilateral solution 

To compute the unilateral solution, we first need to determine whether y U pper (/■*> ^) occurs in the 
iteration for n even or odd. To this end, we temporarily denote the two points as y dd (a^ 8) and 
y e ven (a 2 ?^)- in the unilateral case, i.e. y od \, yivln > 0, equation <JT7|) leads to 

2/odd = A yeven + (/X + (5) C, y e ven = A y odd + (fJL - 8) C. (19) 

Solving the above equations yields 

2/odd (/■*, ^) = £f P (A 4 ) + ^ 2/even (m> S) = x fp (p) - 5d, (20) 

where Xf p and e? are computed in ((TT]) and (|18|). respectively. By adopting the sign convention 

sign (5 = sign (i 1 - 1 ), (21) 

we can arrange that y up per and yi ower correspond to odd- and even-numbered beats, respectively; 
that is 

2/u P pcr (/i, 5) = Xfp (At) +$d, yi owcv (p, 5) = Xfp (ax) - 5 d. (22) 
In the remainder of the paper, we retain the notation y U pp er and yiower i assuming the sign convention 

mm. 
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The unilateral solution is valid provided J/jjjL. > 0, that is 



fi> p\5\= fi CTit (5) , (23) 

where 

> 0. (24) 



y-(l) 



At the critical value of /i, we have 



2/upper (Merit (5) , S) = Xf p (/i crit (<$)) + 5 d = y u _ cr it (S) (25) 
Mower (Merit (*) , S) = x {p (/x crit (<5)) - 5 d = y\_ crit (5) . (26) 



Note that y^rit > and y£ it = 0. 
3.2 The bilateral solution 

We now compute the bilateral solution that evolves continuously from the unilateral solution 1 as 
fj, crosses /i C rit (5). Thus, for \x < //crit (S), we look for a period-two orbit of (|17j) such that 

ySper(^»> >yi ( owe r ^^)' ( 27 ) 

where y U pper occurs for n odd. By Q17|). these satisfy 

2/upper = -B 2/lower + (/X + 5) C, yi ower = A y up per + (M ~ S) C. (28) 

Comparing these equations with l|14|). we see that the solution of ((2*5)) may be written as 

Supper upper Supper ) 2/lowcr — 2- lower (m)~ " Slower j 

(29) 

where x U pper and xi OW er have been computed in Q15j). and the //-independent shifts s up per, Slower are 
given by 

supper = (I-BA) - 1 {I- B) c, s\ ower = {I - AB)- 1 {I - A) c. (30) 

We now show that this solution is consistent with the assumption Q27|). We claim that when 
M = Merit {$), the solution (|25"26() of (|19j) for this value of \x also satisfies (|28j) : this claim follows 
from the observation that, since yp^rit = ^> we nave -^2/i-crit = ^Z/l-erit- Using the fact that x U pper 
and xi ower are linear in /x, we may rewrite (|29j) as 

2/uppcr — Supper (Merit (<5) , ^) + •'-upper (m)- upper (Merit (S)) 

= 2/u-crit (#) + (M - Merit (#)) -^upper (31) 

and 

2/lower (M) <?) = 2/l-crit ($) + (M ~ Merit (<5)) A~i ower . (32) 

It now follows from Condition 6 that (|27|) holds for all /i < /i cr it (5). 

1 In the appendix, we show that for fi sufficiently large negative there is also a second bilateral solution. 
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3.3 Stability of perturbed iterations 

The unperturbed map without alternate pacing as in (|17|). may suffer a so-called dangerous 
bifurcation ^Uj at fi = 0. To elaborate: suppose that at some positive value of \x the (unperturbed) 
iteration is locked onto Xf p (/x). If \x is decreased quasistatically but remains positive, the solution 
will follow the stable fixed point Xf p (/i). However, as /x — > 0, the basin of attraction of the fixed 
point might shrink to nil. In this case, the iteration becomes indeterminant as \x crosses 0, i.e., as 
a result of infinitesimal perturbations, the solution might shift to a different attractor for \i < 
rather than to the period-two solution. 

One may hope that iterations of p7|). i.e., the map under alternate pacing, are better behaved 
near /i = /x C rit (^) when they cross the boundary plane ja^ 1 ) = 0}. However, we are unable to 
guarantee this without imposing an additional hypothesis, as follows. 

Condition 8: There exists an invertible matrix S s.t. ||5'j4 2 5 _1 || < 1 and ||5^4S5 _1 || < 1. 
In other words, the eigenvectors for A 2 and AB are not too different. For example, this condition 
is satisfied if both A 2 and A B are both diagonalizable using S. 

In the following lemma, we consider the composition of (|17|) with itself when \x ~ /i cr it (5), 
starting with n even from an x near the even iterates yi OW er of the period-two orbit, for which 
we have yj^Lr O- 4 ' ^) ~ ^l-crit = Thus, both options in (|TTj) must be considered for the first 
application of ([170. However, the image of x will be close to the odd iterates, for which we have 
2/upper (/•*, S) ~ 2/u-crit > 0- Thus, provided is sufficiently close to fi cr it (5) and x is sufficiently close 
to y^gj. (^,S), only the upper option in (fT7|) will occurs in the second application of (|17jl. Thus, if 
/U and x are so restricted, the composition of ([1711 with itself is given by 

f A[Ax + (/x + 5)c] + (/i-5)c ifx«>0 
W \ A[Ba + 0i + 5)c] + (/i-<5)c ifx«<0 ' 1 j 

Lemma 3.1. F maps a neighborhood of yi ower (/j,,5) into itself and, with respect to an appropriate 
norm, is a contraction there. 

Proof: Given a nonsingular matrix S, we define the 5-norm of a vector v as 

= \\Sv\\ (34) 

and that of a matrix M as 

||M|| S = H5M5" 1 !! . (35) 

For the matrix S in Condition 8, we have = max |||A 2 || S , < 1. 

To prove F in l[33[) is a contraction, we consider x and y in a, neighborhood, of yiower 

and treat 

the following three cases separately. 

1. ) x^\ yW > 0. It follows from (H that 

\\F (x) - F (y)\\ s = \\A 2 (x - y)\\ s < \\A 2 \\ S \\(x - y)\\ s < 9 \\x - y\\ s (36) 

2. ) x( l \ yW < 0. It follows from (j32J) that 

11^ (x) - F (y)\\ s = \\AB (x - y)\\ s < \\AB\\ S \\(x - y)\\ s <9\\x- y\\ s (37) 
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3.) > and < 0. Let z be the point where the line between x and y intersects the 
plane {z^- 1 ' = 0}. This choice of z implies that Az = B z and ||x — z\\ s + ||z — y||g = ||x — y\\ s - It 
follows from (|33|) that 

ll^(^)-^(y)ll s = ll^ 2 *-^yL = ||^ 2 (^-^) + ^5 (*-y)|| s ( 38 ) 

< ||^4 2 || s ||z - z\\ s + \\AB\\ S \\y - z\\ s (39) 
<9\\x- z\\ s + \\y - z\\ s = 9\\x - y\\ s . (40) 

Moreover, since 2/i OW er is a fixed point, F maps a neighborhood of yi OW er into itself.l 

Theorem 3.2. Condition 8 implies that the period-two orbit y U pper, Uiower * s stable. 

Proof: For /x > /i cr ;t and p < /i cr iti stability of the iterations is guaranteed because Vi, 
|Aj (^4 2 )| < 1 and |Aj(l?74)| < 1, respectively. For /i /J c rit; stability is guaranteed because F 
is a contraction, as shown in the proof of Lemma 3.1. ■ 

4 Prebifurcation gain 
4.1 Analysis 

Near the onset of a classical period-doubling bifurcation, disturbances in the bifurcation parameter 
may cause amplified disturbances in the system response, a phenomenon known as prebifurcation 
amplification. Various authors studied prebifurcation amplification for a classical period-doubling 
bifurcation, using a gain, which is the ratio of the response magnitude to the applied perturbation 
amplitude [1611171 ITS] . An analytical formula of the gain, including crucial higher-order terms, was 
derived in ^j]- It is interesting to compute the gain for a B/C period-doubling bifurcation and 
compare with that of a classical period-doubling bifurcation. To this end, let us define a gain using 
the 1st component of the solution 

Y = 2/uppcr Slower 

2 \5\ 

Then, the prebifurcation gain (when p > 0) is described by the following theorem. 
Theorem 4.1. When \i > and 5^0, the gain satisfies 

f r ' const , iffI>p\S\ 

[ r cmS i + 7^-|ffj ifp,<p\8\ ' 

where T con st an d 7 ar & positive constant, and p is given by \24\j - 

Proof: When p > p\5\, the response to alternate pacing is an unilateral solution, as shown in 
the previous section. It follows from Equation Q22|) that 

Ay = 2/upper - 2/lower = 2 5 d. (43) 
By the definition (|41|) and recalling the sign convention (|21|). the gain is 

r= dP> =r const . (44) 
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Figure 3: Schematic of the variation of the canonical gain (|41j) under changes in parameters: (a) T 
vs. p when 5 stays constant and (b) T vs. 5 when p stays constant. 

Therefore, if p is varied and 5 is held fixed, the gain stays constant as long as p > p cr i t (5) = p\5\, 
see Figure 01 (a). Similarly, if 6 is varied and p is held fixed, the gain stays constant as long as 
H < S CTit (p), where 5 CTit (p) = p/p, see Figure El (b). 

When p < p \S\, the response to alternate pacing is a bilateral solution, as shown in the previous 
section. It follows from Equations (|31|) and (|.32|) that 

Ay 

— Uwpper Slower 

= J/u-crit — J/l-crit + (^crit (5) ~ ^) (Supper ~ Slower) (45) 
= 2 5 d + (/i cr it (<5) — ^i) (^upper — Slower) , 

where we have used (|25|) and (|26|) for y u -crit arid j/i_ cr it - Recalling the definition fill) and using (fl4"|) . 
the gain is 

r = r const + 7 (p-^Y (46) 

where 

7 = ^ (4pU - *£L) • (47) 

Here, 7 > 0, since X U p pcr > and X^ cr < 0. Therefore, for a constant (5, the gain increases as p 
decreases from p cr i t = p\S\, see Figure|3](a). On the other hand, for a constant p, the gain increases 
as \S\ increases from |(5 cr it| = p/p, see Figure |3] (b).l 

The precise relations (|42jl . which hold for all p and 5, result from the normal form in 
which higher-order terms have been dropped. If such terms are present, then ()42l) will describe 
the qualitative behavior of the gain near the bifurcation point for small perturbation amplitude. 
In Section 4.2 we present a numerical example with the B/C period-doubling bifurcation in the 
cardiac model of Sun et al. JHj; this example shows that the leading-order analysis accurately 
capture the behavior of the model. 

Remark: Analogous to T, a generalized gain can be defined using an arbitrary component of 
V as 

« _ W 

Ji) _ l/upper l/lower / Aa \ 

9 = 2|*j • ( } 
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It follows from (g3J) and (@SJ that 




gSnst, if/i>p|<5| 

£ ) ns t+(p-fi)^ ) ifM<pi*r 1 j 

where 

<£Lt = (signed) d« fc« = i (xgU " 4L) • (50) 

As we saw above, fi^onst = r con st > and A^ 1 ) = 7 > 0. However, if i ^ 1, the signs of fl^nst 
and feW are not uniquely determined. Thus, the generalized gain may have either sign and, even 
if positive, may exhibit monotonicity opposite to that in Figure 3. However, all gains are constant 
when jj, > p \5\. 



4.2 A numerical example 

Sun et al. |15j presented a theoretical model for atrioventricular nodal conduction, which accurately 
predicts a variety of experimentally observed cardiac rhythms. It was shown in jllj that the model 
of Sun et al. exhibits a border-collision period-doubling bifurcation. Following the notation of 
the model iterates the atrial His interval, A, and the drift in the nodal conduction time, R, as 
follows 

if A n < 130 

(51) 

if A n > 130 

where A m \ n = 33 ms, r rec = 70 ms, Tf a t = 30 000 ms, and 7 = 0.3 ms. The bifurcation parameter, 
H, represents the time interval between bundle of His activation and the subsequent activation. 
Under variation of H, a border-collision period-doubling bifurcation occurs at iJbif = 56.9078 ms, 
where of course A = 130 at the bifurcation point. We apply alternate pacing to the above model by 
perturbing H with (— l) k 5. Figure |1] shows the relation between the gain T and the parameters H 
and 5, computed from our theory as well as from numerical simulations. The agreement between 
the theoretical and the numerical results is good. 



Rn+l 



( ( A min + R n+1 + (201 - 0.7 A n ) e- H ' T - 



A mhi + R n+ i + (500 - 3.0A n ) e-^ec 
R n e~ {An+H),Tt ^ + 7e _H / Tfat 



5 Summary and discussion 

A border-collision (B/C) bifurcation occurs when the fixed point of an iterated map encounters a 
border on which the underlying map is continuous but its derivatives jump. Using the normal form 
(^Q), we presented the conditions for a B/C period-doubling bifurcation in multiple dimensional 
maps. Under these conditions, we computed the response of the map under a subharmonic per- 
turbation to the bifurcation parameter, i.e. alternate pacing. The response under alternate pacing 
can be either unilateral or bilateral. As the names suggest, the two points of a unilateral solution 
lie on the same side of the border, while those of a bilateral solution lie on different sides of the 
border. The response to alternate pacing is quantified through the definition of a gain. We found 
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(b) h - Hut = 
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Figure 4: Alternate pacing of the model of Sun et al. (|51[): (a) variation of the gain T under changes 
of H when (5 = 1 ms and (b) variation of the gain T under changes of 5 when H — = 0, 2, and 5 
ms. Here, the solid curves indicate theoretical results and the dots indicate numerical simulations. 



f-H 









Figure 5: Schematic of the variation of the gain of a classical bifurcation under changes in param- 
eters: (a) T vs. (j, when S stays constant and (b) T vs. 5 when \x stays constant. 

that the gain is a piecewise smooth function of the bifurcation parameter (//) and the perturbation 
amplitude (|<5|), with qualitatively different behaviors in different parameter regions. 

Most importantly, the gain in a border-collision bifurcation differs qualitatively from that in 
a classical bifurcation. For example, it follows from that, after appropriate rescaling of the 
parameters, the gain of a classical period-doubling bifurcation satisfies the relation, 

5 2 r 3 + ^r-i = o. (52) 

Figure (a) and (b) schematically show the dependence of T on /i and T on 6, respectively. 
Comparison of (|42|) and (|52[) reveals two differences between the gain in a classical bifurcation 
and that in a border-collision bifurcation. (1) The gain of a classical bifurcation tends to infinity 2 
as (/j,, 5) approaches (0,0); by contrast, the gain of a border-collision bifurcation is bound for all 
5^0. (2) The gain of a classical bifurcation varies smoothly under changes in system parameters 
while that of a border-collision bifurcation undergoes a nonsmooth variation as parameters cross a 
boundary in the parameter space. 

2 Moreover, the rate of divergence as the parameters (/z, 5) tend to (0, 0) depends on the path taken. For example, 
when 5 is extremely small, the gain tends to infinity as on the other hand, when fj, = 0, the gain tends to infinity 
as5" 2 / 3 . 



11 



Although the two cases exhibit these differences, if T is plotted as a function of [i (with 5 held 
fixed), it may be difficult to distinguish them based on data from experiments, because only discrete 
points can be sampled and these are subject to experimental errors (cf. Figures (a) and 0(a)). 
However, if T is plotted as a function of 5 (with fi held fixed), the distinction between the two 
bifurcation types is evident even with just a few data points and in the presence of experimental 
noise (cf. Figures El (b) and 0(b)). Indeed, the theory in this paper guided the experiments in 
identifying the type of a bifurcation in a paced cardiac tissue. 
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Appendix: The out-of-phase bilateral solution 

Recall that, for the unilateral solution under alternate pacing, y upP er occurs for n odd due to the 
sign convention 1)21(1. and this behavior continues for the bilateral solution computed in Section 3.2; 
i.e. this bilateral solution is in phase with the unilateral solution. Thus, we refer to such bilateral 
solution as an in-phase bilateral solution. On the other hand, when \i is sufficiently large negative, 
there exists a bilateral solution with the opposite phase; i.e., y U pper occurs for n even. We refer to 
the latter as an out-of-phase bilateral solution. 

Continuing the notation as in (|27j). the out-of-phase bilateral solution is characterized by the 
following equations 

2/upper = B yi ower + (p - 6) c, yi owcr = A y npper + (fl + 5) c. (53) 

Note the signs of the terms ±5c are different in the above equations and l)28J) . Solving the above 
equations yields 

2/upper upper 

0*)- 

Supper > J/lowcr Slower Slower ) 

(54) 

where Supper and xi ower have been computed in Equation ()15j) . and the //-independent shifts s U pper, 
si OW er have been computed in Equation l)30|) . Note that the out-of-phase and in-phase bilateral solu- 
tions are parallel to each other because they are both parallel to the solution x uppe r ((a), x\ oweT ((a). 
One can show that the out-of-phase bilateral solution exists only when (aX^ ppct ^ $ Supper- 

As 

required in Condition 6, A upper < 0. Therefore, (i has to be sufficiently large negative for the 
out-of-phase bilateral solution to exist. 
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